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1.  INTRODUCTION  AND  LITERATURE  SURVEY 
1.1  INTRODUCTORY  COMMENTS  ON  THE  UPSTREAM  TRAVELING  WAVE 

In  a  spatial  boundary-value  problem  with  vorticity  introduced  along  a 
grid  in  the  freestream,  Rogler  and  Reshotko  (Ref.  1)  found  an  explosively 
growing  solution  that  represents  vorticitv  diffusing  upstream.  The  solution 
is  not  an  eigenmode.  It  oscillates  neutrally  in  the  direction  perpendicular 
to  the  freestream  rather  than  vanishing  far-away. 

The  counterpart  to  this  solution  appears  in  the  spatial  boundary— value 
analyses  of  Refs.  2-7  (for  2-D  fluctuations)  and  Ref.  8  (for  3-D  fluctuations) 
as  a  branch  line.  The  corresponding  Four ier-Laplace  solution  for  the  velocity 
of  form  0(y)exp[iox-ib>t]  has  a  complex  wavenumber,  u=a _+iaj  (with  a^  very 
large  and  negative,  indicating  rapid  growth),  and  a  real  frequency,  <■>.  It 
represents  a  third  class  of  growing  solutions  of  the  Orr-Sommerf eld  equation. 
This  solution  appears  to  have  no  counterpart  in  the  viscous  temporal 
initial-value  formulations  (Refs.  9  and  10)  with  a  real  and  u>  complex. 
Neither  does  it  appear  in  the  inviacid  spatial  solution  of  Ref.  11. 

The  other  two  classes  of  growing,  spatial  solutions  include  the 
Tollmien-Schl ijhting  (T-S)  wave,  which  is  the  fundamental  eigenmode  that  can 
be  unstable  or  stable  depending  on  the  parameters  and  boundary  layer. 
Although  there  are  other  eigenmodes.  Mack  (Ref.  12)  demonstrated  that  the 
higher  temporal  eigenmodes  are  finite  in  number  and  heavily  damped  for  the 
Blasius  layer. 

Another  growing  solution  is  a  standing  wave.  It  does  not  travel,  but 
grows  exponentially  in  the  stzeamwise  direction  and  oscillates  in  time  as 
expI+fJx-iwt]  with  p  real.  The  phase  speed  is  pure  imaginary.  The  standing 
wave  has  been  described  in  Refs.  13  and  14;  a  movie  was  made  to  visualise  its 
behavior.  The  broader  picture  of  the  five  spatial  classes  of  solutions  is 
presented  in  the  next  section. 

For  high  Reynolds  numbers,  the  upstream  traveling  wave  is  very  difficult 
to  calculate  because  of  extremely  high  frequency  oscillations  in  the 
solutions.  The  high  frequency  oscillation  for  the  upstream  traveling  wave  has 
a  y-wavenumber  O(R^),  as  compared  with  0(Rj  )  for  the  T-S  wave.  At 
*5  =  2500,  this  wavenumber  is  fifty  times  greater  than  for  the  T-S  wave. 

Also,  these  oscillations  can  decay  very  slowly  away  from  the  wall,  with  a 
decay  rate  0(u)  rather  than  0(Rj  )  for  the  instability  wave.  That  means  that 
the  high  frequency  oscillation  can  survive  far  from  the  wall  for  the  upstream 
traveling  wave.  In  contrast,  the  viscous  solution  for  the  T-S  wave  is 
significant  only  near  the  wall  (where  it  represents  the  visoous  sublayer)  and 
the  critical  layer. 

For  these  reasons,  the  present  calculations  have  been  restricted  to  low 
Reynolds  numbers.  At  high  Reynolds  numbers,  this  solution  will  be  a  challenge 
for  future  developments  in  numerical  solutions  of  differential  equations. 

Our  objectives  here  are  to  describe  solutions  of  the  Orr-Sommerf eld 
equation  for  upstream  traveling  vorticity.  Three  cases  are  described: 
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(1)  Analytical  solutions  for  vorticity  in  a  uniform 
flow  far  from  any  wall  (Sections  1.3, 1.4) 

(2)  Analytical  solutions  for  vorticity  in  a  uniform 
flow  near  a  wall  (Section  1.5) 

(3)  Numerical  solutions  of  the  Orr-Sommerf eld 
equation  for  boundary  layers  (Chapters  2  and  3) 

We  believe  that  closely  related  waves  traveling  upstream  also  appear  in 
free-shear  layers,  jets,  wakes  and  other  shear  layers. 


1.2  WHY  STUDY  THE  FOURIER- LAPLACE  SOLUTIONS  OF  THE  ORR-SOMMERFELD  EQUATION? 

To  accurately  simulate  boundary  layer  transition  in  test  facilities,  the 
processes  leading  to  transition  must  be  known  and  controlled.  To  understand, 
cai~”*at.,  measure  and  control  boundary  layer  transition  and  mixing  in  shear 
layeis,  it  is  important  to  identify  all  possible  Four ier-Laplace  solutions  of 
the  viscous,  parallel-flow,  linear  equations.  More  sophisticated  equations 
exist,  but  the  importance  of  these  equations  cannot  be  denied.  They  can 
describe  many  basic  phenomena  in  shear  layers.  To  have  a  mathematically 
complete  set  of  solutions  which  could  be  superimposed  to  satisfy  appropriate 
upstream  and  downstream  boundary  conditions  would  be  a  practically  useful 
accomplishment.  Perhaps  most  importantly,  even  after  more  than  70  years  since 
the  original  derivations  by  Orr  and  Sommerfeld,  and  after  more  than  50  years 
since  the  first  solution  by  Tollmien,  basic  features  continue  to  be  discovered 
from  the  equations  for  viscous,  linear,  parallel-flows. 

Apart  from  the  usefulness  of  a  complete  set  of  solutions  in  initial-value 
problems,  it  is  also  important  to  be  aware  of  physical  phenomena,  to  be  able 
to  categorize  them,  and  to  be  able  to  describe  them  physically.  The 
Fourier-Laplace  solutions  help  sort  out  these  basic  features. 

While  reasons  for  seeking  the  Fourier-Laplace  solutions  are  summarized 
above,  a  few  statements  concerning  what  these  solutions  are  NOT  are  also 
warranted.  In  many  cases,  a  superposition  of  members  of  one  family  of  waves 
or  several  families  of  waves  is  required.  The  behavior  observed  can  be  unlike 
that  of  any  one  Fourier  component,  even  though  that  Fourier  component  is 
imbedded  within  the  observed  behavior.  Secondly,  nonparallelism  and 
nonlinearity  can  introduce  other  significant  features  into  the  problems.  Yet 
even  for  these  cases,  the  Fourier-Laplace  solutions  based  on  the  viscous, 
parallel-flow  equations  provide  (1)  the  solutions  which  can  be  superimposed 
and  (2)  solutions  serving  as  a  point  of  departure  for  weak  nonparallelism  and 
weak  nonlinearity. 

Section  1  of  Ref.  3  surveys  the  temporal  and  spatial  initial-value  and 
boundary- value  problems.  Section  1.6  of  Ref.  14  lists  the  known 
Fourier-Laplace  solutions  of  the  Orr-Sommerf eld  equation  and  the  other 
viscous,  parallel-flow,  linear  equations  for  incompressible  flow. 

The  five  classes  of  2-0  spatial  solutions,  including  decaying  ones,  are 
tabulated  in  Table  I  below.  For  simplicity,  the  table  has  been  restricted  to 
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2-D  fluctuations  (with  u  and  v  velocities,  but  no  w  velocity)  and  a  flat 
plate.  However,  the  pattern  of  solutions  appears  to  generally  repeat  itself 
for  3-D  fluctuations,  i.e.  the  poles  and  branch  lines  from  2-D  formulations 
appear  to  have  generalized  counterparts  for  3-D  studies.  These 
generalizations  include  longitudinal  vortices  and  at  least  three  other 
structures. 

Table  I  has  three  columns:  for  decaying,  neutral,  and  growing  solutions 
in  the  streamwise  direction.  Row  1  is  composed  of  the  set  of  eigensolutions 
which  vanish  far  away  from  the  boundary  layer.  The  fundamental  instability 
wave  (the  Tollmien-Schlichting  wave)  can  be  unstable.  As  indicated  by  the 
arrows,  it  also  can  be  neutral  or  stable  depending  on  the  parameters  and 
boundary  layer. 

Row  1.2  is  a  pair  of  solutions  which  are  vortical  fluctuations  in  the 
freestream.  However,  1.2a  consists  of  vortical  disturbances  propagating 
downstream.  They  are  summarized  in  Section  1.6  of  Ref.  14  and  will  not  be 
surveyed  here.  Although  the  details  of  the  solutions  depend  on  the  problem 
considered,  these  vortical  solutions  appear  in  uniform  flows,  shearing  flows, 
inviscid  and  viscous  flows,  and  in  both  spatial  and  temporal  analyses. 

Solution  1.2b  consists  of  vortical  disturbances  propagating  upstream  and 
growing  very  rapidly  in  the  streamwise  direction.  Fluctuations  1.2b  are 
described  in  this  report. 

The  last  row  also  consists  of  a  pair  of  solutions,  one  which  decays  and 
another  which  grows.  These  two  solutions  are  standing  waves  (they  do  not 
travel  in  the  streamwise  direction)  and  are  irrotational  fluctuations  far  away 
from  the  boundary  layer.  These  two  standing  wave  solutions  are  described  in 
Refs.  13  and  14. 


Table  I 

SPATIAL.  FOURIER/ LAPLACE  SOLUTIONS  FOR  A  BOUNDARY  LAYER 
of  form  fct (y)exp ( iox  -  iwt ) 

for  viscous,  incompressible,  constant  transport  property,  parallel-flow 
over  a  stationary  flat  plate 


Decaying  in  the 
downstream  direction 


Neutral  in  the 
streamwise  direction 


Growing  in  the 
downstream  direction 


(1.1b)  Higher  discrete  modes 

(1.2a)  Downstream  propagating 
freestream  vortices 


(1.3a)  Decaying  standing  waves 


(1.1a)  Tollmien-Schlichting 


(1.2b)  Upstream  propagating 
vorticity  (described 
in  this  report) 

(1,3b)  Growing  standing  wave 
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1.3  2-D  WIVES  DIFFUSING  UPSTREAM  IN  THE  FREESTREAM 

The  Rogler-Reshotko  analysis  of  Ref.  1  is  an  boundary-value  problem  for 
vorticlty  fluctuations  in  the  freestream.  Since  this  analysis  illustrates 
four  classes  of  solutions  in  Table  I  (but  not  the  eigennodes)  with  elementary 
mathematics,  this  analysis  is  suatmarized  in  this  section. 

The  two-dimensional  equation  for  vorticity  fluctuations  in  a  uniform  flow 
is 


?t  +  «x  -  8<$xx  +  5yy)  =  0  (1.1) 

where  e  is  the  inverse  Reynolds  number.  A  generalized  Fourier  transform  in 
the  y-direction  and  time  yields  the  following  equation  describing  the 
evolution  of  vorticity  in  the  streamwise  direction 


filial  (to 


37 


(1.2) 


dx  "  * 

Seeking  a  solution  of  form  £(x)=exp(mx) ,  then  the  characteristic  equation  is 


em1  -  m  +  ( iw  -  ep1)?  =  0  (1.3) 

with  solutions 

■x,s  =  (1  T  [1  -  4e ( iu  -  sp1)]1/*}/2s  (1.4) 

Hence  the  complementary  solution  is 

S 

£(x.3,<i>)  -  C1exp(m1x)  +  Caexp(max)  (1.5) 

The  first  solution  represents  fluctuating  vorticity  which  convects  and  decays 
in  the  streamwise  direction.  The  second  solution  grows  exponentially  in  the 
streamwise  direction.  Both  solutions  can  be  physically  meaningful. 

If  an  initial-value  problem  were  solved  for  the  region  upstream  of  the 
y-axis,  then  the  first  solution  would  be  discarded  and  the  second  solution 
retained.  For  the  region  downstream  of  the  y-axis  where  the  vorticity  is 
specified,  then  the  first  solution  would  be  retained  and  the  second  solution 
is  discarded.  If  a  rectangular  strip  were  considered,  then  both  solutions 
could  exist  in  the  same  domain. 

The  inviscid  and  viscous  solutions  are  sketched  in  Figure  1.1.  Without 
viscosity,  vorticity  cannot  diffuse  upstream.  The  solution  is  identically 
zero  until  the  fluid  crosses  the  y-axis  where  the  vorticity  is  introduced. 
The  vorticity  is  convected  downstream  without  decay. 

With  viscosity,  however,  vorticity  can  diffuse  in  the  upstream  direction. 
Instead  of  the  (inviscid)  Heaviside  distribution,  two  exponential  functions 
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n, — 

envelope  of  fluctuating  vorticity 

upstream 
decay  length 
vanishes 

X 

downstream  characteristic  decay  length  is  infinite 

1.1a  Inviscid  case  with  a  uniform  mean  flow 


1.1b  Viscous  case  with  a  uniform  mean  flow 


1.1c  Effects  of  reducing  the  mean  velocity  while  holding  the 
dimensional  frequency  and  y-  and  z-wavenumbers  constant. 


Figure  1.1  Envelopes  for  the  inviscid  and  viscous  cases  of  fluctuating 
vorticity  input  along  the  y-axis. 
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join  at  tlie  y-axis.  The  decay  in  the  upstream  direction  it  "very  rapid", 
since  the  vorticity  mast  diffuse  upstream  against  the  mean  flow.  In  the 
downstream  direction,  the  decay  is  "slow"  as  the  vortices  are  convected 
downstream  and  diffuse. 

Both  viscous  solutions  may  be  relevant  in  both  regions,  if  both  positive 
and  negative  mean  flows  arise  in  both  regions.  An  example  is  the  case  of 
reversed  flows.  Another  example  is  the  region  between  the  leading  and 
trailing  edges  of  a  finite-length  plate. 

Generalizations  of  the  second  solution  to  account  for  a  boundary  layer 
axe  an  objective  of  this  research. 


1 .4  3-D  WAVES  DIFFUSING  UPSTREAM  IN  A  UNIFORM  FLOW 

The  equation  describing  the  evolution  of  small-amplitude  disturbances  in 
a  uniform  flow  (U=l)  reduces  to 


l*  +  t  -  iv*KA  =  0  (1.6) 

dt  ox  Eg 

Letting  v=0(y)exp [iox+iYz-i»t] , 

then  the  Orr-Sommerf eld  equation  for  a  uniform  mean  flow  reduces  to 

U  -  “  -  JL_(D*  -  a2  -  t*)](D*  -  a*  -  r*)l*  -  0  (1.7) 

a  iaRg 

The  characteristic  length  6  in  the  nondimensionalizations  is  a  dummy  parameter 
for  this  uniform  flow  case.  While  that  length  could  be  defined  as  the 
characteristic  viscous  length,  the  same  nondimensionalizations  will  be  used 
here  as  defined  in  Chapter  2  so  that  our  analytical  results  can  be  compared 
with  the  numerical  solutions. 

For  solutions  which  oscillate  sinusoidally  in  the  y  direction 

0(y)  *  CcosPy  +  DsinfJy  (1.8) 

the  Orr-Soaimerf eld  equation  reduces  to 

U  -  “  +  *  (a*  +  p2  +  r*)](o*  +  p*  +  r*)  =  0  (1.9) 

a  iaRg 

This  equation  has  four  solutions  for  the  x-wavenumber,  a.  Two  solutions  are 
irrotational  standing  waves  which  grow  and  decay  in  the  streamwise  direction 
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a  =  ±i(p*  +  y*)1/a  (1.10) 

Th«se  wave*  are  described  in  Ref.  13  and  14. 

Two  other  solutions  are  rotational,  traveling  waves 

a  =  ~i*§.{l  -  [1  +  _i_(p*  +  r*  -  iwR6)]l/*)  (1.11a) 

V 

a  =  "i*6{l  +  [1  +  _J_(p*  +  y*  -  iuRg)]1/*}  (1.11b) 

*  R8* 

Solution  (1.11a)  with  the  minus  sign  corresponds  to  vortical  fluctuations 
decaying  in  the  streamwise  direction.  These  waves  correspond  to  those  of  1.2a 
in  Table  I.  Solution  (1.11b)  with  the  positive  sign  is  the  one  of  interest  in 
this  report. 

For  large  values  of  the  Reynolds  number,  the  series  representation  of  the 
square  root  can  be  found.  The  results  for  the  x-wavenumber  are 


or  =  -<u [1  -  2e*(p*+Y1+«*) 

+  2e4(p*+r*4«*)(3p*+3y#+7«»*)  +  0(e‘)l  (1.12a) 

«  -i  -  e(0,+7,-We1 ) 

e 

+  e’(p*-ty*4«*)(p*4y*+5i**)  +  0(e‘)  (1.12b) 

where  e  «=  1/Rg.  Based  on  the  first  term  of  the  expansion  for  oj,  the  wave 
grows  in  the  streamwise  direction  as 


exp(-aix)  **  exp(Rgx)  (1.13) 

To  lowest  order,  this  growth  rate  is  independent  of  the  wavenumbers  and 
frequency.  The  solution  for  the  fluctuations  can  be  written  as 


v  =  0<y)  e_aix  eiar<*“cept)eiT* 
where  the  equi-phase  speed,  a  real  quantity,  is 


(1.14) 


cep  =  <a/ttr  =  -1  -  2e*  (p*+y*+»* ) 

+  284(p14T*+«*)(p,4Y*+5»*)  +  0(e‘)  (1.15) 

This  quantity  is  the  speed  of  a  point  of  constant  phase  on  the  wave,  such  as  a 
local  maximum  or  minimum,  a  point  of  zero  value,  etc.  This  speed  is  negative, 
indicating  that  the  waves  propagate  upstream. 
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Table  II  tabulates  a  few  values  of  wavenumber  and  equi-t^ase  speed  for 
various  Reynolds  numbers.  In  the  limit  as  R^  *  <*>,  the  waves  propagate  at  a 
speed  equal  to  the  freestream  speed,  Uw,  in  the  upstream  direction. 


Table  II 


x-Wavenumber s  and  Phase 

Speeds  for  Upstream 
=  w  =  0.5,  y=0 

Traveling  Waves 

*8 

“r 

°i 

cep=«/ar 

1.0 

-0.3217971 

-1.2768869 

-1.553774 

10.0 

-0.4951209 

-10.049272 

-1.0098543 

100.0 

-0.4999500 

-100.00500 

-1 .0000999 

1000.0 

-0.4999995 

-1000.0005 

-1 .0000010 

00 

-0.5 

“*6 

-1.0 

1 .5  WAVES  PROPAGATING  UPSTREAM  IN  A  UNIFORM  MEAN  FLOW  NEAR  A  WALL 

Some  features  of  waves  propagating  upstream  near  a  wall  can  be  extracted 
from  the  uniform-flow  version  of  the  Orr-Sommerf eld  equation.  The 
x-wavenumber  is  given  by  eqn.  (1.11b).  Seeking  solutions  of  form  exp(my),  the 
four  independent  solutions  are 


0  =  Ae  “y  +  Be"* "W  +  CcosPy  +  DsinfJy 
with  B  =  0,  C  =  1,  and  m  =  (a1  +  y*)1^1 

Introducing  the  impermeability  and  no-slip  conditions,  then  A 
and 


(1.16a) 

(1.16b) 

-1,  D  =  -m/p. 


0(y) 


-e-By  +  cospy  - 


"sinPy 


(1.17) 


This  velocity  fluctuation  is  plotted  in  Figure  1.2.  Note  the  high 
frequency  oscillation  which  decays  slowly  away  from  the  wall.  Also  note  the 
undamped  sinusoidal  oscillation  in  the  freestream. 


In  eqn. (1.16a),  the  velocity  has  been  nondimensional ized  by  setting  the 
constant  C=1 .  In  retrospect,  the  choice  D=1  would  be  a  better  choice  for 
large  Reynolds  numbers.  If  the  right-hand  side  of  eqn. (1.17)  were  multiplied 
by  -p/m,  the  solution  with  this  alternate  nondimensional ization  would  result. 


Expressions  (1.16  and  1.17)  are  for  3-D  disturbances.  To  obtain 
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relations  naefnl  for  comparison  with  the  2-D  calcnlations  in  Chapter  2,  the 
following  formulae  ire.  i££  2^  flnctnations. 

With  a  positive  frequency,  the  valne  of  the  exponent  ■  is 


m  =  -u  =  o>  +  iR^  (for  y  =  0)  (1.18) 

The  term  exp(-my)  has  the  following  approximate  valne  for  large  Reynolds 
numbers 

e_By  =  e-iRsy  "  “y  (for  y  =  0)  (1.19) 

Hence,  the  term  exp (-my)  can  be  a  high  y— wavenumber  wave  which  decays  slowly 
in  the  y-direction.  as  shown  in  Figure  1.2. 

The  longitudinal  velocity 


u  =  f(y)eiox"iwt  (for  y  -  0)  (1.20a) 

can  be  found  from  the  continuity  equation.  ux  +  Vy  *  0, 

f(y)  =  +ie-lny  -  *£sinf)y  -  icospy  (for  y  =  0)  (1.20b) 

a 

This  amplitude  is  plotted  in  Figure  1.3  for  R&bIO.  The  high  frequency 
oscillation  in  f(y)  is  proportionately  much  stronger  than  for  0.  Note, 
however,  that  the  scale  of  the  abscissa  has  been  changed.  To  illustrate  the 
effect  of  changing  Reynolds  number,  the  amplitude  f(y)  is  plotted  in 
Figure  1.8  for  Rg=20.  This  figure  clearly  shows  the  higher  frequencies 
associated  with  the  term  exp(-iR-y).  fhile  the  explicit  solutions  of  this 
chapter  can  be  evaluated  for  much  higher  Reynolds  numbers,  the  numerical 
solutions  of  the  Orr-Sommerf eld  equation  presented  in  the  next  two  chapters 
cannot  be  found  for  large  Reynolds  numbers  because  of  the  high  frequency 
oscillation  exp(-iR^y)  which  dampens  slowly  as  exp(-wy)  if  the  frequency  is 
small.  For  this  reason,  all  plots  and  the  numerical  calculations  have  been 
restricted  to  low  Reynolds  numbers. 

From  the  definition  of  vorticity. 


?  =  vx  ~  Uy  =  Z(y)eiax-iwt  (for  y  =  0)  (1.21a) 

the  amplitude  is 
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i  high  frequency  oscillation  can  be  teen 
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Z(y)  =  ~(a  ^  )cos0y  +  a(a  +f ' ) slnpv  (for  y  =  0)  (1.21b) 

la  lap 

This  amplitude  is  plotted  in  Figure  1.4.  Note  that  there  are  uo  high 
frequency  oscillations;  the  term  exp (-my)  is  irrotational.  Also  note  that 
there  is  no  evidence  of  large  amplitude  spikes  of  vorticity  at  the  wall 
associated  with  viscous  sublayers.  The  sinusoidal  oscillations  are  undamped 
and  represent  one  part  of  the  vortical  freestream  disturbances. 

The  fluctuating  pressure 


p  =  Jt(y)eiax  “  (for  y  =  0)  (1.22) 

can  be  found  from  the  x-m omentum  equation 

ut  +  ux  m  ~Px  +  ^(uxx+Uyy)  (for  y  =  0)  (1.23) 

Substituting  for  u,  then 

"(y>  '  ‘I  -  1  -  4lt<yl  *  <f“  >  ‘ 0>  (1-24> 

The  pressure  is  plotted  in  Figure  l.S.  The  presence  of  the  wall  strongly 
influences  the  pressure  through  the  term  exp(-my). 

The  streamlines  and  iso-vorticity  contours  are  plotted  in  Figures  1.6  and 
1.7.  The  influence  of  the  wall,  through  the  high  frequency  term  exp(iKgy)  can 
be  clearly  seen  in  the  streamline  pattern.  A  more  subtle  influence  of  the 
boundary  conditions  at  the  wall  is  that  they  dictate  the  value  of  the  constant 
D  in  eqn. (1.16a).  This  influence  of  the  wall  does  not  vanish  far  away  from 
the  wall.  Only  this  effect  on  the  phasing  of  the  oscillation  remains  in  the 
iso-vorticity  contours  of  Figure  1.7. 
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Figure  1.6  do  not  appear  in  the  equi-vorticity  contours  when  the 
flow  is  uniform. 
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2.  FORMULATION  OF  THE  MATHEMATICAL  SYSTEM  FOR  3-D  FLUCTUATIONS 
INTERACTING  WITH  BOUNDARY  LAYERS 

The  linearized,  parallel-flow,  viscons  momentum  equations  can  be  combined 
into  the  following  equation  in  terms  of  the  fluctuating  velocity,  v(x,y,z,t), 
in  the  direction  normal  to  the  plate 


(t  +  U(y)£.  -  =  0  (2.1) 

art  5x  k6  57  ji 

where  the  Laplacian  is  Va=3a /9xa+9a/3ya+3a /9za .  The  Reynolds  number  based 
upon  the  freestream  velocity  and  characteristic  thickness  of  the  boundary 
layer  is  Rft  =  pU^S/p.  For  Blasius  boundary  layers,  the  characteristic 
thickness  is 

6  -  (2px/ PUW) a^a  (2.2) 

The  lengths  in  eqn.  2.1  are  nondimensionalized  against  6  and  the  time  is 
nondimensionalized  against  i/H^,  We  now  seek  solutions  of  form 

v(x,y,z,t)  =  0(y)exp[iax+iyz-iwtl  (2.3) 

where  the  wavenumbers  and  the  frequency  are  related  by  eqn.  (1.11).  The 
x-wavenumber  is  complex,  while  the  z-wavenumber  and  frequency  are  real. 
Substituting  eqn.  (2.3)  into  (2.1),  then  after  some  rearrangement,  the 
Orr-Sommerfeld  equation  results 


HU  -  “)(Da  -  oa  -  r1)  -  U_  -  1_  (D1  -  oa  -  ya)a)0  -  0  (2.4) 

a  yy  iaRg 

where  D  =  d/dy.  This  equation  is  subject  to  impermeability  and  no-slip  at  the 
wall  and  boundedness  far  away 


O  =  D0  =  O>ty  =  O  and  0  is  bounded  is  y-»s  (2.5) 

When  U  =  1  outside  the  boundary  layer,  the  system  reduces  to 


[1  -  “  -  _^__(Da  -  a*  -  ya)](D*  -  o*  -  r*>J0  -  0  (2.6) 

a  iaKg 

Seeking  solutions  of  form  exp(ssy),  then  the  four  solutions  are 


exp(±ipy)  or  cos£y  and  sinPy  (which  oscillate  neutrally) 
and  exp(lmy)  (which  oscillate  and  grow  or  decay)  (2.7) 

where  m  is  the  root  with  positive  real  part 
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*  -  (a*  +  y%)x/*  (2.8) 

Hence  the  eolation  ontside  of  the  boundary  layer  it 

0(y)  =  Ae~®y  +  Be+my  +  Ccos0y  +  Dsinpy  (for  y  >  y#)  (2.9) 

where  B  =  0  in  the  present  problem  and  we  normalize  by  setting  C  =  1.  Ve 

require  that  the  exterior  eolation  and  the  first  three  derivatives  agree  with 

the  interior  solution  and  the  corresponding  derivatives  at  the  "edge"  of  the 

boundary  layer,  y»=y. 

6 


D^y")  =  D^Xy*)  for  n  «  0,1,2 ,3 


(2.10) 


With  the  normal  velocity  found  as  a  solution  of  the  Orr-Sommerf eld 
equation,  then  the  longitudinal  velocity  is  found  from  the  continuity 
equation.  The  vorticity  is  calculated  from  the  derivatives  of  the  velocities, 
and  the  pressure  is  found  from  the  i -momentum  equation,  rms  quantities  and 
other  correlations  such  as  the  Reynolds  stress  and  energy  production  are  also 
found. 


Our  objective  now  is  to  solve  the  mathematical  system  composed  of  the 
Orr-Sommerfeld  equation  (2.6a,b),  the  two  wall  conditions  (2.8),  and  the  four 
conditions  at  the  boundary  layer  edge  (2.10). 


25 


AEDC-TR-83-7 


3 .  NUMERICAL  SOLUTIONS  FOR  WAVES  TRAVELING  UPSTREAM  IN  BOUNDARY  LAYERS 


3.1  SUMMARY  OF  THE  NUMERICAL  TECHNIQUES 

The  numerical  techniques  for  calculating  both  the  mean  flow  and  the 
unsteady  flow  are  the  sane  as  used  in  Refs.  13  and  14.  The  numerical 
techniques  are  summarized  in  this  section. 

To  obtain  the  coefficients  U(y)  and  U  (y)  of  the  mean  flow,  numerical 
solutions  of  the  Falkner-Skan  equation  'are  obtained  by  the  method  of 
Nachtsheim  and  Swigert  (Ref.  15.  The  equation  is  solved  by  shooting  from  the 
wall  outward,  using  4th-order  Runge-Kutta  integration  in  double-precision 
arithmetic.  The  exterior  boundary  condition,  F  •  1  as  q  •■»,  is  not  specified 
at  the  boundary  layer  edge,  y  -  ye,  but  rather  the  error  is  minimized  in  a 
least  square  sense. 

The  Orr-Sommerfeld  system  of  equations  (2.5b,  2.8,  2.13,  2.14)  is  solved 
by  an  expansion  in  Chebyshev  polynomials 

0((|)  -  (3.1) 

SH*0 

with  polynomials  defined  as 

TB(q)  *  costm(arccosq)]  (3.2) 

and  where  the  variable  q  is  related  to  the  dimensionless  y  by  the  linear 
transformation 

n  -  -  1  (3.3) 

A  set  of  linear,  algebraic  equations  results  from  the  Orr-Sommerfeld  equation, 
the  two  wall  boundary  conditions,  and  the  four  matching  conditions  at  the 
boundary  layer  edge.  The  matrix  of  coefficients  for  this  system  is  reduced  by 
Gauss-Jordan  elimination  with  maximum  pivoting  by  column  and  rows.  With  the 
coefficients  known,  the  expansion  3.1  is  carried  out  for  0(y),  and  identities 
used  to  obtain  the  y-derivatives  of  0(y) .  With  0  and  its  derivatives  known, 
the  longitudinal  velocity,  vorticity,  and  vorticity  production  are  calculated. 
From  the  x -moment urn  equation,  the  pressure  fluctuations  are  found,  and  then 
the  correlations  are  calculated  for  the  rms  quantities.  The  streamlines  and 
iso-vorticity  contours  are  also  drawn. 

The  numerical  checks  include  numerical  experiments  to  determine  the 
effects  of  number  of  terms  in  the  expansion,  precision  of  arithmetic  for  the 
matrix  reduction,  and  influence  of  the  y-value  where  the  edge  boundary 
conditions  are  imposed.  Comparisons  are  made  between  analytical  and  numerical 
solutions  for  a  uniform  mean  flow.  The  same  computer  program  used  to  obtain 
solutions  for  the  traveling  waves,  except  with  different  patching  conditions 
at  the  boundary  layer  edge,  was  also  used  to  recover  the  Tollmien-Schl ichting 
waves,  the  two  standing  waves,  decaying  vortical  fluctuations  in  the 
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freestream,  and  decaying  longitudinal  vorticea. 

3 .2  NUMERICAL  RESULTS  FOR  UPSTREAM  TRAVELING  WAVES  IN  A  BLASIUS  BOUNDARY  LAYER 

The  velocity  profile  and  derivatives  for  the  Blasint  boundary  layer  are 
plotted  in  Fignre  3.1  with  the  same  y-scale  as  the  other  figures. 

To  aid  in  the  interpretation  of  these  calculations,  since 


v  =  [0r(y)+i01(y)]elarx-i“te-aix 


=  0rcos(arx-e*t)e-aix  -  0£sin(arx-«t)e’’aix  (3.4) 

then  for  urx-«t  =0,  v  =  0r(y))e-aix 
«rx-«t  =  n/2,  v  =  -0i(y))e-aix 
«rx-«t  =  n,  v  =  -0t(y) )e-aix 

<*rx-<i»t  =  3n/2,  v  =  0£(y))e-aix  (3.5a-d) 

The  vertical  lines  urx-vit=0 ,  jt/2,  it,  and  3n/2  travel  to  the  left  at  speed  cep 
=  w/ur  (slightly  less  than  -1) .  Bence  the  real  and  imaginary  parts  of  the 
complex  amplitude  (when  multiplied  by  the  factor  expt-a^x))  are  the  velocity 
profiles  along  those  vertical  lines  arx-«t  =  0  and  3n/2  respectively.  With  a 
reversal  of  sign,  they  are  the  velocity  profiles  along  the  lines  arx-wt  =  n/2 
and  n.  Analogous  interpretations  also  apply  for  the  other  complex  amplitudes. 

Figures  3.2  through  3.15  are  plots  for  a  wave  with  y-wavenumber  and 
frequency,  p=u=  0.5,  and  Reynolds  number  Rg«  10.  These  calculations  are  for 
2-D  fluctuations,  although  the  computer  program  for  0(y)  is  applicable  for  3-D 
fluctuations. 

In  Figure  3.2  for  the  normal  velocity,  note  that  the  impermeability 
condition  is  satisfied  at  the  wall.  Also  note  the  dual  oscillatory  behavior 
of  the  solution.  Far  from  the  boundary 


0(y)  =  cosPy  +  Dsinpy  (3.6) 

This  figure  (with  a  Blasius  layer)  can  be  compared  with  Figure  1.2  (with  a 
uniform  mean  flow).  Note  the  substantial  influence  of  the  boundary  layer  on 
both  the  magnitude  of  the  high-frequency  oscillations  and  the  phasing  of  the 
low-frequency  oscillations  (sin0y  and  cosf)y)  which  survive  into  the 
freestream. 

Figure  3.3  is  a  plot  of  the  rms  normal  velocity  fluctuation  which  shows 
that  impermeability  is  satisfied  and  illustrates  the  slow  dampening  of  the 
high-frequency  oscillation  in  the  y-direction. 


In  Figure  3.4  for  the  longitudinal  velocity,  note  that  the  no-slip 
condition  is  satisfied  and  again  note  the  general  oscillatory  behavior  of  the 
solution.  This  figure  (with  a  Blasius  layer)  can  be  compared  with  Figure  1.3 
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Figure  3.4  Calculation  of  the  longitudinal  velocity  fluctuation  for 
P  =  si  =  0.5,  y  =  0,  R&  =  10  and  a  Blaaiua  boundary  layer.  Thia 
figure  (for  a  Blaaiua  layer)  can  be  compared  with  Figure  1.3  (for 
a  uniform  mean  flow). 
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(for  a  uni for*  mean  flow).  While  both  figures  illustrate  the  dual  oscillatory 
behavior,  the  details  of  the  flows  are  quite  different.  The  y-wavenumber, 
frequency,  and  Reynolds  number  are  the  same  for  both  figures. 

Figure  3.5  is  the  rms  version  of  the  previous  figure  which  illustrates 
the  slow  dampening  of  the  high-frequency  oscillation. 

The  streamlines  for  the  fluctuating  flow  at  time  zero  are  plotted  in 
Figure  3.6.  Again  note  the  dual  oscillatory  patterns  in  the  y-direction  and 
the  effects  of  the  exponential  growth  in  the  streamwise  direction,  behaving 
approximately  as  exp(+R6x).  When  comparing  Figure  3.6  (with  a  Blasius  layer) 
and  Figure  1.6  (with  a  uniform  mean  flow),  note  the  much  stronger  high 
frequency  oscillations  in  Figure  3.6. 

The  fluctuating  vorticity.  5=v  -Uy,  is  plotted  in  Figure  3.7;  the  rms 
vorticity  is  plotted  in  Figure  3.8.  In  contrast  to  Figure  1.4  where  no  high 
frequency  oscillations  are  present.  Figures  3.7  and  3.8  clearly  show  that  such 
oscillations  arise  when  the  boundary  layer  is  present. 

The  source  of  the  high  frequency  rotational  fluctuations  when  the 
boundary  layer  is  present  is  shown  in  Figure  3.9,  which  is  a  plot  of  the 
vorticity  production  term  in  the  equation  describing  the  evolution  of  vortical 
fluctuations  in  the  boundary  layer 

?t  +  °5x  +  vly  =  *<5Xx  +  tyy)  (3.7) 
where  v?y  =  -  0(y)Uyy(y)exp(iax-i«t) 

This  term  is  analogous  to  the  term  describing  the  ^production  of  disturbance 
energy  in  shear  layers,  -uvdU/dy.  The  ter*  -f>0yy  is  plotted  as  Figure  3.9 
which  shows  that  the  high-frequency  velocity  produces  high  frequency 
vorticity.  The  vorticity  which  has  been  produced  appears  in  Figures  3.7  and 
3.8  as  the  high  frequency  vortical  fluctuations  near  the  wall.  This 
production  vanishes  at  the  wall  (since  v=0  there)  and  in  the  freestream  (since 
0yy=0  outside  of  the  boundary  layer). 

The  iso-vorticity  contours  are  plotted  in  Figure  3.10.  This  pattern  is 
substantially  simplier  than  the  streamline  pattern  of  Figure  3.6,  but  contains 
the  high-frequency  oscillation  which  Figure  1.7  does  not  have  when  the  mean 
flow  is  uniform. 

The  complex  amplitudes  and  the  rms  amplitude  of  the  fluctuating  pressure 
are  plotted  in  Figures  3.11  and  3.12.  Much  of  the  high-frequency 
characteristics  of  the  instantaneous  pressure  is  lost  by  the  rms  version,  at 
least  for  this  particular  case.  Note  that  the  pressure  is  not  constant  across 
the  boundary  layer. 

The  Reynolds  stress,  -uv ,  and  the  kinetic  energy  production,  -uvdU/dy, 
are  plotted  in  Figures  3.13  and  3.14.  The  Reynolds  stress  again  reflects  the 
slowly  decaying,  high-frequency  oscillation.  The  energy  production  is 
confined  to  the  boundary  layer.  It  vanishes  at  the  wall  where  v  =  0  and  in 
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Figure  3.5  Calculation  of  the  ms  longitudinal  velocity  fluctuation  for 
8=u>  =  0.5,y=0,  R.  *=  10  and  a  Blasius  boundary  layer. 


Calculation  of  the  vorticity  fluctuation  for  P  =  u  =  0.5 
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the  freestreaa  where  Dy  =  0. 

The  averaged  kinetic  energy  (u  +v  )/2  of  the  nneteady  flow  is  plotted  in 
Figure  3.15,  which  illustrates  that  a  region  exists  in  the  boundary  layer 
where  the  kinetic  energy  is  larger  than  in  the  freestreaa.  This  region 
roughly  corresponds  to  the  region  where  the  energy  production  is  positive  in 
Figure  3.14. 
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4.  SUMMARY,  DISCUSSION  AND  CONCLUSIONS 
4.1  INTRODUCTORY  COMMENTS 

Thi*  report  describes  a  fluctuation  which  propagates  upstream  in  a 
boundary  layer.  This  fluctuation  diffuses  in  the  upstream  direction,  againat 
the  mean  flow,  and  dampens  very  rapidly  in  the  upatream  direction.  The 
propagation  speed  is  always  more  negative  than  -U„  (i.e.  the  absolute  value 
of  the  speed  is  always  greater  than  UM  and  is  directed  upstream) .  Aa  the 
Reynolds  number  increases,  the  The  speed  approaches  -Uw  as  the  Reynolds  number 
becomes  infinite. 

These  fluctuations  oscillate  sinusoidally  in  time.  Par  from  the  boundary 
layer,  they  oscillate  sinusoidally  in  the  y-direction. 

Evidence  of  the  existence  of  these  waves  is  based  on  numerical  solutions 
of  the  Ort- Sommer f eld  equation,  analytical  solutions  of  simplified  forms  of 
the  Orr-Sommerfeld  equation,  and  an  asymptotic  solution  of  the  Orr-Sommerfeld 
equation.  The  asymptotic  solution  is  not  soomarized  in  this  report,  but  a 
similiar  solution  was  used  by  Shunichi  Tenge'  (Ref.  Id)  for  the 
Tollmienr-Schlichting  wave.  The  use  of  that  asymptotic  solution  in  the  spatial 
initial-value  problem  is  sonmarized  in  Refs.  2-4. 

It  is  instructive  to  outline  three  equations: 

(1)  The  viscous  equation 

It +  "  s^1)v‘v  ‘  V*>Jr  - 0  (4*1) 

has  a  4th-order  derivative  in  x 

4th-order  derivative  in  y 
lst-order  derivative  in  t 


This  equation  yields  four  basic  solutions  far  from  the  boundary 
layer:  The  decaying  and  growing  vortical  fluctuations  (which 

propagate  downstream  and  upstream  respectively)  and  the  decaying  and 
growing  standing  waves  which  are  irrotational .  The  family  of 
eigenmodes  vanish  far  from  the  wall. 

(2)  The  inviscid  form  of  the  above  equation 

l*  +  U(y)*  )y2v  -  U  (y)*I  =  0  (4.2) 

at  3x  3x 

has  a  3rd-order  derivative  in  x 

2nd-order  derivative  in  y 
lst-order  derivative  in  t 
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Far  fro*  the  boundary  layer,  this  equation  yields  only  three  basic 
solutions:  A  nondecaying  traveling  wave  of  vortlcity  (representing 
vortices  being  convected  downstream),  and  a  pair  of  decaying  and 
growing  standing  waves.  The  solution  corresponding  to  the  present 
work  is  missing. 

(3)  The  viscous  equation  for  the  case  of  no  mean  flow  reduces  to 


{£.  -  1v*KA  =  0  (4.3) 

«Ft  U6 

and  has  the  same  order  of  derivatives  as  case  (1) ,  with  three 
solutions  corresponding  to  vortlcity  diffusing  downstream  and 
upstream,  and  two  standing  waves.  In  this  case,  the  exponential 
decay  of  the  two  viscous  solutions  are  symmetric:  the  asymmetry 
introduced  by  the  mean  flow  and  x-dependence  through  the  term  Ud/dx 
is  lost.  A  similar  case  arises  for  soanwiae  fluctuations  for  the 
special  case  of  vanishing  x-wavenumber . 


4.2  SUMMARY  OF  THE  FORMULATION  AND  NUMERICAL  TECHNIQUE 

For  small-amplitude  disturbances  in  shear  layers  which  develop 
sufficiently  slowly  such  that  the  parallel-flow  representation  is  adequate, 
the  objective  is  to  obtain  solutions  of  the  equation 

it  +  my)t  -  Iva)V*v  -  D  (y)il  *  0  (4.4) 

<Ft  <Tx  Rg  yy  dx 

which  are  of  the  Fourier-Laplace  form 

v(x,y,z,t)  =  0(y)exp[iax+iyz-i»t]  (4.5) 

where  the  frequency  u  is  real,  and  the  x-wavenumber  is  complex 

1  +  [1  +  ^.(p*  +  r*  -  io*6)]l/*}  (4.6) 

The  motivations  for  seeking  numerical  solutions  of  this  form  are  that  such 
waves  appear  in  analytical  solutions  where  the  mean  flow  is  uniform,  as  shown 
in  Chapter  1.  Furthermore,  in  an  asymptotic  analysis  of  an  initial-value 
problem  with  a  boundary  layer,  a  branch  line  indicated  that  waves  of  this  form 
would  exist  in  boundary  layers.  For  solutions  of  this  form,  the  equi-phase 
speed  is 
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c#p  =  »/ar  =  -1  -  2s* (0*+1T *■*»*) 

+  284(p,+r*+«*)(p*+T*+5«1)  +  0(e‘)  (4.7) 

The  impermeability  and  no-slip  boundary  conditions  at  the  nail,  and  the 
boundedness  condition  of  the  solution  far-away,  are 


0*D0**Oaty  =  O  and  0  is  bounded  as  y  ■»  »  (4.8a,b,c) 


When  U  -  1,  the  Orr-Somaerf eld  equation  has  the  solutions 


0(y)  ■  Ae-By  +  Be+®y  +  Ccosfly  +  Dsinpy  (for  y>ye)  (4.9a) 

or  equivalently 


0(y)  =  Ae  *y  +  Be+my  +  Caexp(-i0y)  +  Daexp(if)y)  (for  y  >  ye)  (4.9b) 

where  the  exponent  is  a  *  (u*  +  y*)1^*  (4.10) 


For  a  boundary  layer,  B  =  0.  For  aost  cases,  solution  (4.9a)  is  noraalized  by 
setting  0=1 . 

Far-awav  froa  the  boundary  layer  edge  where  the  term  Aexp(-ay)  is 
negligible,  the  flowfield  varies  sinusoidally 


0(y)  “  Ccosfiy  +  DsinfSy  (for  y  »  ym)  (4.1D 

This  £xi9ttX11M  disturbance  &  rotational.  The  constant  D  is  obtained 
nuaerically. 

In  general,  the  constants  A  and  D  and  the  solution  0(y)  of  the 
Orr-Soaaerf eld  equation  depend  on  the  frequency,  y-wavenuaber,  Reynolds 
nuaber,  and  the  mean  velocity  profile. 

Numerical  solutions  of  the  Orr-Sosuaerf eld  equation  at  low  Reynolds 
numbers  are  obtained  by  an  expansion  in  a  series  of  (typically  48)  Chebyshev 
polynomials.  The  coefficients  of  those  polynomials  are  obtained  as  solutions 
of  a  set  of  linear,  algebraic  equations.  The  astrix  of  coefficients  is 
reduced  by  Gauss-Jordan  elimination.  Calculations  are  carried  out  on  a 
Digital  Equipment  Corporation  PDP-11/23  computer  using  64K  bytes  of  memory  and 
another  64K  bytes  of  extended  memory  for  virtual  arrays.  The  disk-overlayed 
FORTRAN  IV  program  executes  in  about  145  seconds  when  using  the  KEF11-AA 
floating  point  chip,  beginning  with  data  for  the  mean  velocity  profile,  and 
ending  with  the  coefficients  of  the  polynomials  and  the  amplitude  0(y) .  The 
prograa  executes  in  about  80  seconds  when  accessing  the  FPF-11  floating  point 
board. 
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4.3  SUMMARY  OF  RESULTS,  DISCUSSION  AND  CONCLUSIONS 

Hie  analyses  and  calculations  indicate  that  a  very  rapidly  growing 
solution  of  the  Orr-Sommerfeld  equation  exists  which  differs  in  form  and 
behavior  from  the  Tollmien-Schlichting  wave  (the  fundamental  eigensolution) . 
The  wave  has  the  usual  Fourier-Laplace  form,  eqn. (4.5).  The  oscillation  grows 
very  rapidly  as  exp(+Rgx)  in  the  streamwise  direction.  The  solution  is  a 
traveling  wave  which  propagates  upstream  with  a  speed  greater  than  the 
freestream  speed,  U^,  i.e.  cep<-l.  In  the  inviscid  limit,  the  growth  rate 
becomes  infinite.  This  result  reflects  the  idea  that  diffusion  upstream 
cannot  occur  without  viscosity.  The  rapid  decay  in  the  upstream  direction 
arises  because  the  vorticity  must  diffuse  upstream  against  the  mean  flow. 
This  growing,  upstream-traveling  wave  also  differs  in  form  from  the  growing 
standing  wave  of  Refs.  13  and  14. 

In  the  uniform  mean  flow  outside  of  the  boundary  layer,  the  solutions 
oscillate  as  a  superposition  of  neutral  sinusoidal  waves,  cos0y  and  sin0y,  as 
well  as  another  high-frequency  oscillation.  Because  this  high-frequency 
oscillation  may  decay  very  slowly,  it  can  survive  into  the  freestream.  In  the 
freestream,  this  high-frequency  oscillation  is  irrotational.  Inside  the 
boundary  layer,  high-frequency  rotational  fluctuations  are  generated. 

Since  the  solutions  admit  oscillations  of  the  form  exp(-iRgy),  the 
calculations  have  been  restricted  to  low  Reynolds  numbers  where  evidence 
exists  of  convergence.  By  increasing  the  number  of  terms  in  the  series 
expansion,  and  perhaps  also  the  precision  of  the  calculation,  higher  Reynolds 
number  cases  could  be  calculated.  The  author  notes  that  the  computer  program 
which  has  calculated  Tollmien-Schlichting  waves  at  RgMOOO  has  only  marginal 
convergence  at  Rg=10  for  the  present  wave  with  48  terms  in  the  series 
expansion.  For  this  reason,  the  numerical  solutions  plotted  in  Chapter  3  are 
preliminary. 

In  Refs.  13  and  14,  the  role  of  the  growing  standing  wave  in  numerical 
solutions  of  the  partial  differential  equations  was  discussed.  That 
discussion  is  repeated  below,  since  the  growing,  rotational  solutions 
documented  in  this  report  will  also  arise  in  numerical  solutions. 

Refs.  17-24  are  examples  of  numerical  solutions  which  have  been  obtained 
for  unsteady  boundary  layer  problems  with  finite  difference,  vortex  filament, 
and  spectral  techniques.  These  investigators  have  been  faced  with  the 
difficulty  of  properly  posing  the  boundary  conditions  on  the  downstream, 
outer,  and  upstream  boundaries  of  the  calculational  regions.  This  problem 
with  boundary  conditions  will  become  more  severe  as  the  disturbances  adopt  the 
characteristics  of  freestream  turbulence  and  its  interaction  with  the  leading 
edge,  as  the  amplitudes  increase,  and  as  3-D  and  other  significant  features 
are  included. 

The  growing  vortical  wave  could  be  a  factor  in  obtaining  numerical 
solutions  in  a  rectangular  domain,  since  those  waves  are  a  mechanism  by  which 
the  downstream  boundary  condition  can  influence  the  flow  elliptically. 
Improperly  formulated  boundary  conditions  can  excite  extraneous  fluctuations 
in  the  calculational  domain.  Numerical  errors  from  roundoff  and  truncation 
could  grow  extremely  rapidly  if  a  marching  scheme  in  the  downstream  direction 


.N 

-.'W  ■%."  « 


.vV'-VO 

vv.,-*. 


■V-V.-iV: 

‘‘.'.V.N 

•\«v 


7* 

*,*  *.*  *  v  *•  -  -  h  -  -  - .  *,  •.  -\  -V  /.  >.  V  V  V  v  . vN  -  * .  v  v\*  .  •  ■*  ..■*  .v  .• j 


AEDC-TR  83  7 


were  employed  and  if  this  solution  were  excited  and  uncontrolled.  El initiation 
of  these  dangerous  waves  may  be  necessary. 

In  other  cases  as  suggested  by  the  paragraphs  above,  the  proper  boundary 
conditions  should  incorporate  vortical  waves  which  naturally  enter  into  the 
calculational  region  across  the  upstream  boundary.  This  wave  is  a  viscous 
mechanism  for  upstreaa  communication  in  this  elliptic  problem. 

Analogous  waves  also  appear  in  3-D  cases.  For  at  least  one  of  the  cases 
of  spanwise  fluctuations  in  a  boundary  layer,  the  analogous  diffusive  wave 
decays  away  from  a  vorticity  source  much  less  rapidly  than  the  streaawise  case 
documented  here. 
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NOMENCLATURE 


coefficient  of  the  mth  Chebyshev  polynomial 

constants  in  the  uniform-flow  solution  of  the 
Orr-Soaimerf eld  equation 

complex  phase  speed 

equi-phase  speed 

ordinary  derivative  in  the  y-direction 
nth  ordinary  derivative  (D°<>=©) 

averaged  disturbance  kinetic  energy  per  unit  mass 
complex  amplitude  of  the  longitudinal  velocity 
(-l)l/l 

exponent;  index  on  Chebyshev  polynomial 
number  of  terms  in  the  series  expansion 
disturbance  pressure 

Reynolds  number  based  on  characteristic  thickness  of  the 
boundary  layer 

time 

mth  Chebyshev  polynomial 

disturbance  velocities  in  the  x,y  and  z  directions 
mean  x-velocity  in  the  streamwise  direction 
mean  velocity  in  the  freestream 

coordinate  parallel  to  plate  and  in  streamwise  direction 
cordinate  normal  to  the  plate 
y-value  of  the  boundary  layer  "edge" 
lateral  or  spanwise  coordinate 


Z(y)=Zr+iZi 


complex  amplitude  of  the  fluctuating  vorticity 
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Greek  end  Script 
a 

P 

r 

v* 

6 

e=l/R6 

n 

c > 

5 

n(y)*=nr+ini 

0(y)=Or+i®i 


x-wavenumber 
y-wavenumbe  r 
z- wave number 

8*/8x1+8l/3y*+81/8z*  3-D  Laplacian  operator 
characteristic  thickness  of  the  boundary  layer 
inverse  of  Reynolds  number 
(2y/ye)-l 

kinematic  viscosity 

disturbance  vorticity  in  the  z-direction 

complex  amplitude  of  pressure  disturbance 

complex  amplitude  of  the  normal  velocity  disturbance 

streamf unction 

frequency 


Superscripts,  Subscripts,  and  Miscellaneous  Notation 
(  )  the  time  average  over  one  time  period 

e  boundary  layer  edge 

r.i  real  and  imaginary  parts  of  a  complex  variable 

<>x  partial  derivative  with  respect  to  x 

A  double  generalized  Fourier  transform 
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